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§1. Introduction 

Elliptic regularity theory is traditionally concerned with functions / : 
Q — > IR, £1 being some domain in some Euclidean space IR m , or, more gener- 
ally, in a Riemannian manifold M. f is assumed to solve some elliptic PDE, 
and the regularity theory derives estimates of various norms of / in terms of 
some data (boundary conditions, geometry of the domain) and some struc- 
tural constants of the elliptic operator. The prototype here is the Laplace 
operator, and the elliptic equation in the simplest case then reduces to the 
Poisson equation 

Af = v, 

for some prescribed v, plus some boundary condition. The idea of elliptic 
regularity theory then is to control some Sobolev or Holder norm of /, 

\\f\\ W k, Pl Or ||/|| C fe,a 

in terms of the corresponding norm of A/ for k — 2 in place of k plus some 
terms depending on the data. This may then be applied to estimate a 
solution of the Poisson equation through the prescribed right hand side v. 
As is well-known, such estimates provide the basis for the existence theory for 
solutions of elliptic equations and a guide for the construction of numerical 
approximation schemes (for a detailed presentation, see e.g. [3]). 

We are interested here not in functions, but in maps / : Q — ► N into some 
Riemannian manifold N. This will make the problem genuinely nonlinear. 
Nevertheless, the Laplace operator A/ naturally generalizes to the tension 
field r(/), given in local coordinates by 
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with the standard notation: 

(7a/3( a; ))a,/3=i,...,m is the Riemannian metric tensor on the domain !lcM 
w.r.t. local coordinates x = (x 1 , . . . , x m ), (7 a ^)oi,/3=i,..., m its inverse, 7 = 
det(7 Q( g). r*- fe (/) are the Christoffel symbols of the Riemannian metric tensor 
{9ij{f))i j=l,...n on N, where we take the liberty to identify / with its local 
coordinate representation (/, . . . , f n ). 

Those local coordinates can be chosen in such a manner that the Christof- 
fel symbols r^ fe (/) vanish precisely if N is flat, and in that case, the tension 
field reduces to the Laplacian of the domain f2. In general, however, a Rie- 
mannian manifold is not flat, and so the tension field then is a nonlinear 
elliptic operator. 

In more abstract terms, r(/) is a section of the bundle f~ l TN, and 
therefore fundamentally nonlinear as this bundle itself depends on the map 

This tension field is the negative gradient field of the energy functional 

E{f) = \J \df(x)\ 2 dvol (M) = y^( x) g tJ (f( x ))^W^dx 1 ...dx 

in the same way that the Laplace operator is the negative gradient field of 
the Dirichlet integral. Critical points of the energy integral, i.e. solutions of 

r(f) = 0, 

are called harmonic maps and have been intensively studied. 

More generally, given some function : N — > IR, one may look for critical 
points of the integral 

E*(J) :=\J \df(x)\ 2 dvol (M) + J *(/(*)) dvol (M). 

Solutions solve a system of type 

r(f)=V(f), 

for some vector field V on N, and are called harmonic maps with potential, 
see e.g. P, 0, Ell- 
in the present paper, we shall be concerned with the system 

r(f(x)) = V(x,f(x)), 
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for some given V, without assuming a variational structure, however. Since 
r(/) is a section of f _1 TN, also the right hand side naturally has to depend 
on /, in contrast to the linear Poisson equation where V is a function of x 
only. 

The existence problem has been studied by von Wahl in case of a 
domain 0, with boundary and prescribed boundary values, and by Jost-Yau 
[6] for this case as well as for the more subtle one of a closed manifold M. 
In particular, it turns out that, without a variational structure, on a closed 
manifold M, the above problem need not possess a solution. This is in 
contrast to the Dirichlet problem on 0, that has a solution under general 
circumstances, as shown in those papers. 

Actually, all these results need some curvature assumptions on the target 
N, and, without imposing a size restriction, one needs the assumption that 
TV" have nonpositive sectional curvature. Therefore, we shall assume through- 
out this paper that -/V has nonpositive sectional curvature. Under this 
assumption, we wish to study the above problem as a nonlinear generaliza- 
tion of the Poisson equation. Still, while we no longer have a linear structure, 
the curvature assumption implies that there exists an underlying convex ge- 
ometry as has been explored in more abstract terms in [I]. 

The first part of our paper is concerned with extending elliptic regularity 
theory to the present nonlinear setting. The guiding idea is that the tension 
field r should assume the role of the Laplacian A. An additional feature, 
however, is that not only the geometry is nonlinear, but also that the topol- 
ogy is nontrivial. Therefore, our estimates will also involve a topological 
datum, namely the homotopy class of the map in question. With that addi- 
tion, we are able to generalize the fundamental elliptic W 2,2 estimate to our 
setting. In Theorem. 2.2, we shall estimate the L 2 -norms of the first and 
second derivatives of any map / : M — > N solely in terms of the L 2 -norm of 
its tension field plus a topological term that only depends on the homotopy 
class [/] of /. 

Of course, the constants in those estimates will depend on the underlying 
geometry, in particular on a bound for the Ricci curvature of M. 

In the second part of our paper, we study the Dirichlet problem for 

r(f) = V. 

Our result will need some bound on V depending on the first eigenvalue of 
ft. The result as well as the method are different from those obtained by 
von Wahl and Jost-Yau. 
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§2. A nonlinear W 2,2 — estimate 



A continuous map / : M :— > N induces a homeomorphism 

P = f# : 7Tl(M) — > 7Tl(JV) 

of fundamental groups and at the same time a lift 

f : M —> N. 

We shall need the p— equivariance of the lift map /, i.e. 

/(As) = p(A)/(x) 

for all x G M, A G 7Ti(M) where the fundamental groups m (M) and 7Ti (A/") 
act isometrically on M and TV by deck transformations respectively so that 
M = M /ir\{M) and N = N/tti(N). There is a correspondence between the 
p— equivariant maps from M to N and maps from M to N . 

For the complete simply connected Riemannian manifold of nonpositive 
sectional curvature N, the distance function 

d : N x N — > ]R, u ) = the distance from u to t> 

is well defined and we have a smooth function d 2 on the manifold N x N. 
Let /i, /2 : M — > TV be p— equivariant maps, then the function d 2 (fi,f2) is 
also smooth on the manifold M because deck transformations are isometric. 
For (u, v) G TV x N, we choose an orthonormal basis &\, &2, ■ ■ ■ , G n for T U N. 
By parallel translation along the shortest geodesic from u to v on TV, we 
get a basis e±, e2, ■ • • , e n for T^/V. Take ei, e2, ■ ■ ■ , e n , ei, e2, ■ ■ ■ , e n as a local 
orthonomal frame for Ti u ^ v \N x N. Let 6 1 ,9 2 • • ,6 m be an orthonormal 

coframe in a neighbourhood of x G M, we then write the differentials of the 
maps /i, / 2 : M -»• /V as cZ/i = /^e, ® 6» a , and d/ 2 = /^e^ ® 6» a . We use the 
Einstein summation convention. 

Now we define E(fi, /2), the energy of the difference of the two maps /i, 
h by 



1 /" - 2 

/ 2 ) = g J E E (/L - /L) dvol(M) 

iL/r i=l a=l 



where the integral is calculated on a fundamental domain of M in M. Since 
parallel transport is isometric, we obtain a triangle inequality 



E(h,h)-JE{h,f 2 ) <JE(f 1 ,f 2 ), 
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Inparticular, the energy of the difference is symmetric, and for the energy of 
a map /, 



1 r 2 
^/) = 2/EE(/«) dvol{M), 



M 



i=l a=l 



we then have 



E(f x )-JE{f 2 ) <JE{h,f 2 



(2.0) 



Putting X a = fl a ei + f^a^ii we then are able to express the Laplacian 
of d 2 (f 1 ,f 2 ) as 

Ad 2 (/i, h) = (d 2 )x a x a + d{dJl aa + djflj 

where r{fi) 1 = f{ aa and r(/i) 2 = / 2aa are the components of the tension 
fields of the maps f\, / 2 respectively. It was shown in ^U] that if Kn < 
then 



Hence 



Ad 2 (/!,/ 2 )> 2]T (fia-f: 

i 

We have also 



2o 



d{\r(h)\ + \r(f 2 )\}. 



Ad(/i,/ 2 )>-{|r(/i)| + |r(/ 2 )|}. 



(2.1) 



(2.2) 



Let /i and / 2 be maps from a domain of the manifold M, smooth enough in 
Q and on the boundary 90. If two maps fx and / 2 coincide on the boundary 
<9f2, we then have by (2.2) 



\Vd{fuf2)\ 2 dvol < ||rf(/i,/2)|| La(n) {lr(/i) 



Therefore 



L 2 (n) 



+ 



L 2 (Q). 



^^ 2 )|L 2 (n)^ A ^"{lk^lL 2 (n) + I^IL 2 (n)} M 



L 2 (C) 



where the constant A(f2) is the first eigenvalue of the Laplacian on Q, see 
On the other hand, the integral of (2.1) over Q gives 

E(f 1 ,h)<2- 1 j Ad 2 (/ 1 ,/ 2 ) + d(/ 1 ,/ 2 ){|r(/ 1 )| + |r(/ 2 )|}. 
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Notice that eZ 2 (/i,/ 2 ) > in Q and d 2 (fi,f2) = on dQ, the derivative 
of d 2 (/i,/2) along the outer normal vector of £2 at d£l will be nonpositive. 
Hence the Stokes formula gives 



J Ad 2 (/!,/ 2 )<0. 



We then obtain an energy estimate for two maps that have the same bound- 
ary values and belong to the same homotopy class, with (2.3) 



E(f 1 J 2 )<X(nr 1 {\\r(f 1 



L 2 (Q) 



+ 



T(/2 



2 

L 2 (n). 



(2.4) 



For the case of a compact manifold M without boundary, we need a deeper 
analysis. Integrating the inequality (2.1) over a fundamental domain, we get 



E{h,f 2 )< \ d(h,f 2 ) 



L 2 (M) 



{||r(/i 



L 2 (M) 



+ 



r(/ 2 



L 2 (M) 



}. (2.5) 



A direct corollary of the triangle inequality and the inequality (2.2) is 
that if /ii,/i2 are two homotopic harmonic maps, then E(h\, h,2) = and 
E{fM) = E(f,h 2 ) whenever [/] = 

Lemma 2.1: Zei /o, /i : M — > TV" 6e p— equivariant maps. Define 
ft : M ^ N by exponential map 

ft(x) = exp Mx) {texp- ( ^ x) fi(x)}, 



then ft is also p— equivariant and yjE(ft) is a convex function oft. 

Proof: Let ^ x be the geodesic on TV from fo(x) to fi(x), then ^ x is also 
p— equivariant. Hence, ft is p— equivariant. Fix a € {1,2, • • • , n}, we claim 
that V a ft(x) is a Jacobi field along the geodesic^. In fact, let u a G T X M be 
the vector dual to 8 a , c(s) be the geodesic on M with c(0) = x, c'(0) = v a , 
C : (0, 1) x (-6,6) -► TV", C(t,a) = exp /o(c(s)) {t exp^ 1 ^ /i(c(s))}, then for 
any fixed s, C(-, s) is a geodesic. On the other hand 



s=0 
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By the Jacobi field equation, we have 



^^E\U) = -^= J (R(V a f t (x),i) i,V a f t (x))dvol(M) 



+ 



^y/E||v a / t (x)| 2 ^(M) 



i 



M 



E(ft 
> 0. 



Vaft(x)V a ^- t f t (x)dvol(M)\ 2 



M 



As required. 

Lemma 2.2: There exists a constant C which depends only on the ho- 
motopy class [f] such that 



d(f,h) 



L 2 (M) 



<CE(f,h) 



(2.6) 



holds for a harmonic map in the homotopy class [f] . 

Proof We prove the Lemma by contradiction. If there were a homotopy 
class such that the inequality (2.3) did not hold, we then could find a se- 
quence fk of maps and harmonic maps h k in the same homotopy class such 
that 



= inf 

L 2 (M) T (h)=0 



d(fk,h) 



L 2 (M) 



with 



d(fk,h k ) 

2 

l 2 (m; 

is independent of the choice of the harmonic map h. Define 

flix) = exp hk(x) {texp^ (x) f k (x)}. 



>k 2 E(f k ,h k ) (2.7) 



d(f k ,h) 



> 1. Notice that the energy of the difference E(f k , h) 



We have 



2 

L 2 (M) 



d(f k ,h 



d(fl,h k ) 
>k 2 t 2 (^Ejfk)-^E(h k )) 2 . 



L 2 (M) 



>k 2 t 2 E(f k ,h k ) 



It follows from the convexity of the energy E(fj t ) that 
d(f k ,h k ) _ , > k 



L 2 (M) 



E{f k ) - ^E{h k ) 



(2.8) 
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Choose t = tk such that d(fl, h k ) 



L 2 (M) 



1. From (2.8), fj, k is a minimiz- 



ing sequence for the energy. Therefore the sequence converges strongly 
to a harmonic map h, i.e. 



d(f k \h) 



L 2 {M) 



0. 



(2.9) 



Notice that 



d(fkJl k ) 



L 2 (M) 

(1 -t k )d(fk,h k ) 



d(f k ,h k ) - d(h k ,f k k ) 



L 2 (M) 



d(fk, h k ) 
d(fk,h k ) 
Hence for k large enough, 



L 2 {M) 
L 2 (M) 



L 2 (M) 

d(f k , hk] 



- tk 

- 1. 



L 2 (M) 



d(fk,h) 



L 2 (M) 



< 


d(fk 


ft) 




d{fk 


hk) 



L 2 (M) 
L 2 (M) 



1 

+ 2 
1 

~ 2' 



This is a contradiction to the choice of the harmonic maps /ifc. This proves 

2 

> 1. Notice that 



L 2 (M) 



that the inequality (2.3) holds for inf T (^) =0 d(f, h) 

the harmonic map and the inequality (2.3) are invariant under the rescaling 
of the metric on N , the inequality (2.3) holds for all maps in the same ho- 
motopy class. 



Now we summarize the above, that is, (2.0), (2.5) and (2.6) as: 
Theorem 2.1: Let M be a compact Riemanian manifold with or with- 
out boundary, N a compact Riemannian manifold of nonpositive sectional 
curvature. For a given homotopy class of maps from M to N, with prescibed 
boundary values when the boundary of M is not empty, there is a constant 
C such that for any map f in this homotopy class there is a harmonic map 
h in the homotopy class such that 

(j \df?f <(J \dh?f +C(j \r{f)?f (2.10) 

MM M 
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where r(/) is the tension field. 

The first term on the right hand side depends only on the homotopy class 
of /, because all harmonic maps in the same homotopy class have the same 
-minimal- energy. 

Recall the Bochner formula ([H], Ch. 8) 

Ae(/)(x) = \Vdf\ 2 + (8(T(f)),df) + {df(Ric M (e a )),df(e a )) 
- (R N (df(e a ), df{ep))df{e p ),df{e a )). 

Integrating it on the domain M and making use of the inequality (2.10) 
under the assumption of nonpositive curvature, we get the main result of 
this section, namely an estimate for the L 2 — norm of the first and second 
derivatives of a map in a given homotopy class in terms of its tension field. 

Theorem 2.2: Let M be a compact Riemanian manifold without bound- 
ary, N a compact Riemannian manifold of nonpositive sectional curvature. 
For a given homotopy class of maps from M to N, let h be a harmonic map 
in that homotopy class. We then have, for any map f in that homotopy 
class, 

J \df\ 2 dvol{M)+ J \Vdf\ 2 dvol(M) < C x J \T(f)\ 2 dvol(M)+C 2 J \dh\ 2 dvol(M) 

M M M M 

where the constant C\ depends only on the homotopy class of the map f and 
the constant C2 is 1 + H-ffic^Hoo, where Ric M is the Ricci curvature of M. 

We point out once more that the last term on this inequality depends 
only on the homotopy class (and on the geometry of M and N), but not in 
any way on the map / in that homotopy class that we are estimating here. 

By differentiation, we may then also obtain higher order estimates in a 
standard manner. 



§3. Boundary value problems 



We now let fi be a domain in a manifold M with a nonempty boundary 
<9f2 and compact closure 0. For the moment, we assume that the map 
g : n -> N is of class C 2 ' a . 

Consider the parabolic system 



r(f) = V(f), 

f(x,0)=g(x), xeQ, 
f{x 2 t) = g(x), x E (90, 
/ : Q x [0, 00) -> N. 



(3-1) 
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Remark: If there is an underlying variational structure, as for harmonic 
maps with potential, namely, 

EM) = \J \df\ 2 dvol(M) + J <p{f{x))dvol{M), 
n n 

for some function ip : N — > R, we have the Euler-Lagrange equation as 

r(/)(x) = V^(/(x)). 

We have 

^E v (f(;t)) = -J \^f(x,t)\ 2 dvol(M), 
n 

for a solution of the associated parabolic problem 

^=r(/)-Vp(/), 
and therefore, if we assume that E v (f) is bounded from below, i.e., 

> c 

for some C € IR and all /, then subconverges to zero in I? for 

i — > co, and the analysis becomes easy. This has been explained in the 
literature, see e. g., (HI- Without a variational structure, however, the 
problem is more difficult. 

The smoothness of the map g and the theory of linear parabolic systems 
give us a short time solution of the parabolic problem (3.1), i.e., there is a 
positive T so that there is a C 2,a solution on Q x [0, T). The condition we 
impose on the vector field is that 

X7V(X,X) > -fi\X\ 2 , XeT(Tn). (3.2) 

We also let \(Q) be the first nontrivial eignvalue for the Dirichlet problem 
on S7. 

Lemma 3.1: Let f satisfy the parabolic system (3.1) where the vector 
field V satisfies (3.2) with fi < %\(U), then 

S\y{-M<S\r{9)-V{g)\\ (3.3) 
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Proof: A direct computation gives that 

U 

if 



-2(^(M,|{)/^,f). 



(3.4) 



With the assumptions of the Lemma and the nonpositivity of the sectional 
curvature of the manifold N, it follows from the above identity 



Notice that 



( A - W 



V 



di 



> -§A(J2) 



3/2 dj_ 

dt ' dt 



> 



at 



V 



+ 2 



v dt 



(3.5) 



dt 



Multiplying the two sides of the inequality (3.4) by 
over Q, one obtains 



ld_ 

2dt 



9/4 > 3 
dt ~ 2 



l dt 1 



-A(fi) 



and then integrating 



df 



dt 



Because ^ = on the boundary of the domain $7, the right side of the above 
inequality is nonnegative. Hence 



< 0. 





df 


dtj 


dt 



This proves Lemma 3.1. 

Let us return to the inequality (3.5) 



a df d_l c dj_ df 
(A dt K dt>dt } - C{ dt'dt } 



Set ip(x,t) = exp{Ci}(§£,§£), then 



(A-|)^>0. 



(3.6) 



and ip = on the boundary of S7. By Moser iteration (cf |S]), one obtains 
that 

V> 2 (Mi) <c(l + — ?— ) / V 2 (M) 

with the constant C independent of t. It follows from Lemma 3.1 and the 
above inequality that 

df 4 



dt 



<C \ T (g)-V(g)\' 



11 



Therefore we get a uniform bound for 1 1£ | . Now we have 

Ad(f(;t),g(.))>-C (3.7) 

where the constant C comes from the uniform bound of and the vec- 
tor field V by (2.2) and the parabolic system (3.1). A direct corollary of 
the inequality (3.7) is that the distance function d(f(-,t),g(-)) has a uni- 
form bound which is independent of t. Indeed applying the weak maximum 
principle to the elliptic inequality (3.7) one gets 

supd(/(-, t), <?(•)) < sup <*(/(-, t), </(•)) +c\nfi. 
n an 

Consider the Dirichlet problem for the inhomogeneous Laplace equation 
in 0, 

f Au = -C, in Q; . . 

1 u = o, on an, { ' 

where the constant comes from (3.7). 

Let u be the solution of problem (3.8), then 

A{d(f(-,t),g(-))-u(-)}>0. 

Hence 

supH/o, -«(■)} < supM/O, t), <?(•)) -«(■)}■ 

n an 

Therefore 

d(f(-,t),g(-)) <u(-). (3.9) 

The solution of the problem (3.8) is of class C 2 ' Q if the boundary dQ 
smooth enough. So we have 

u(x) < Cd{x,dQ). 

This implies 

d(f(;t),g(-)) <Cd(x,dn). (3.10) 

Take v±, ■ ■ ■ , v m to be a local orthonomal frame at z$ £ <9S1 with i/i, • • • , f m -i 
the tangent vectors of 90 and z^ m the normal direction of the boundary. At 
the point zq G <9f2 we then have 

°JM = m, !<,<„-!, 
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by the boundary condition, and 

df(;t) 



< C 



because of (3.10). Therefore we obtain a uniform bound for the gradient of 
the maps /(•, t), i.e., 

e(f)(x) <C, xedn 
where the constant C is independent of t. 

^From the above argument we know there is a constant R so that 
d(f(-,t),g(-))<R on a 
Let F = 2R 2 — d 2 {f{-, t), £/(•)), then F > R 2 . Set 

We have 

(A 9 )A(r ^- ( A -MW) 2e(/)(A-|)F 4Ve(/) ■ VF 6e(f)\VF\ 2 

(3.11) 

By the Bochner formula, (cf 0, Ch. 8), 

(A-|)e(/) = ((A-|)^-/,^- / ) 

= (V/„„ t v; + (v„ iWJ /, v ViU j) , 312 , 

>\B(f)\ 2 -Ce(f), 

where B(f) = V FiVj / is the Hessian of the maps f(-,t). On the other hand, 

(A - |)F = V 2 N F(f*u h Uui) - (V N F, r(f) - f ) 
= V 2 N F{Uv i ,Uv i )-{V N F,V). 

By Jacobi field estimates, (see, e.g. (0], Ch. 4) 

V%F(f m Vi,Uvi) < -2e(/). 

Hence 

(A-^-)F<-2e(f) + 2R\V\. 
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Here we have used the property that \VnF\ = 2d{f{-,t),g{-)) < 2R. Re- 
turning to (3.11), we have 

(A - §- t )A(x,t) > ±e\f) - ^e(f) - ^e(f) + /, (3.13) 



where 



/ = jr 2 W)! 2 " ^Ve(/) • VF + ^\VF\\ 



Because 

„ , FVe - 2eVF 

VA = J3 ■ 

we have 

/= ></)!"- ^|Ve (/ ).V F+ ^|V^-^. 
Notice that 

\B(f)\ 2 - ||Ve(/) • VF + ^|VF| 2 > (\B(f)\ - \df\\VF\) 2 , 
we have 

> 2VA ■ VF 
~ F ' 

Therefore 

( x d \ a, ,w 4 2/- C ... 4fl|F| ... 2VA-VF 
\ dt) A ^ * F3 6 (/) " i^ e(/) " ^ e(/) F ' 



Or 

(A _| M(l , t)a4 ^_^-^-^^. ,3,4) 

Let 

A(t) = max A(x, t), 

and assume that A(t) = A(xt,t) for some point xt- If xt € dfl for some t, 
we have 

e(f){x,t) < 2R 2 A{x,t) < 2R 2 A{x t ,t) < 2e{f)(x t ,t). 
It follows from (3.11) that 

e(f)(x,t)<2C (3.15) 
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where the constant C is the same as in (3.11). On the other hand if xt € ft, 
we have 



--A\ {xut) >4FA'-CA-^A 
Without loss of generality, we may assume that ^-A|( Xt ,t) > 0- Hence 



That is to say 



Therefore 



4FA 2 (t) - CA{t) - ¥*jp.A(t) < 0. 

F 



A(x,«) < £ ' ^ 



4F F 2 



e(/)(x, t) < CR 2 + R\\V\\, (3.16) 

where ||V|| = max ye jv 

Up to now, we have proved 

Theorem 3.1: Let the vector field V satisfy the condition (3.2) with 
fi < | A(J7) for the first eigenvalue A(fi) of the Laplacian in // the maps 
f(-,t) satisfy the parabolic system (3.1) for < t < T, then there is a 
constant C which is independent of t so that 

\df(-,t)\, \%\<C. (3.17) 

By linearizing and using the theory of linear parabolic systems and the 
implicit function theorem, one gets from (3.17) that (3.1) has a solution for 
all of t. That is 

Corollary 3.1: Let the vector field V satisfy the condition (3.2) with 
fi < §A(S7) for the first eigenvalue A(f2) of the Laplacian in Q, then the 
parabolic system (3.1) has a solution for all time t € [0, oo). Moreover, the 
solution of the parabolic system (3.1) has a uniform C 2 ' a bound. 

The last claim of Corollary 3.1 means that any sequence tk — > oo will 
contain a subsequence t n > — > oo so that /(-,£*/) converges to a map / in C 2 . 
This leads to the existence for the Dirichlet problem 

{ = V M- *o (3-18) 
\ f(x) = g(x), xedQ v ' 

under a somewhat stronger condition than Theorem 3.1. 

Theorem 3.2: Let the vector field V satisfy the condition (3.2) with 
H < |A(J7) — e for the first eigenvalue A(f2) of the Laplacian in Q, then the 
Dirichlet problem (3.18) has a solution. 
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Proof: By the same calculation as in the proof of Lemma 3.1, one gets 



that 



ld_ 

2dt 



r 0/4 3 r 
J dt ~ 2 J 



V 



^1 
dt 1 



(lm) -2s) J 



df 



dt 



The Poincare inequality gives 

d_ 
di 



That is to say 



/ 


df 


4 <-4 £ / 


df 




dt 




dt 



n 





df 


1 


dt 



< C exp{-4et}. 



Hence, for a sequence tk with f(-,tk) convergent in C 2 , 

t fe ^oo dt 

Therefore the limit / of /(•, tk) will solve the problem (3.18). 
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